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ABSTRACT. The stability of the bright solitary wave solution to the perturbed cubic- 
quintic Schrodinger equation is considered. It is shown that in a certain region of parameter 
space these solutions are unstable, with the instability being manifested as a small positive 
eigenvalue. Furthermore, it is shown that in the complimentary region of parameter space 
there are no small unstable eigenvalues. The proof involves a novel calculation of the 
Evans function, which is of interest in its own right. As a consequence of the eigenvalue 
calculation, it is additionally shown that A-bump bright solitary waves bifurcate from the 
primary wave. 


1 Introduction 


The nonlinear cubic-quintic Schrodinger equation (CQNLS) is given by 

iAi = Axx + |A|^A + a|A|^A, 


( 1 . 1 ) 


where A is a complex-valued function of the variables (x, t) G R x R+. When a = 0, 
the equation becomes the focusing cubic nonlinear Schrodinger equation, and is used to 
describe the propagation of the envelope of a light pulse in an optical fiber which has a 
Kerr-type nonlinear refractive index. For short pulses and high input peak pulse power 
the refractive index cannot be described by a Kerr-type nonlinearity, as the index is then 
influenced by higher-order nonlinearities. In materials with high nonlinear coefficients, such 
as semiconductors, semiconductor-doped glasses, and organic polymers, the saturation of 
the nonlinear refractive-index change is no longer neglible at moderately high intensities 
and should be taken into account (|l^). Equation (|1.1D is the correct model to describe the 
propagation of the envelope of a light pulse in dispersive materials with either a saturable 
or higher-order refraction index (|l^, |pT| ). 

Equation (^) cannot really be thought of as a small perturbation of the cubic nonlinear 
Schrodinger equation, as it has been shown that a physically realistic value for the parameter 


a IS a ~ 


0.1 (|15|). It turns out that the most physically interesting behavior occurs when 
the nonlinearity is saturating, so for the rest of this paper it will be assumed that a < 0 
(Q, 1^, [1^1, 1^1, [^]). An optical fiber which satisfies this condition can be constructed, 
for example, by doping with two appropriate materials (@). 

One of the more physically interesting phenomena associated with the double-doped 
optical fiber is the existence of bright solitary wave solutions (|A(x)| —> 0 as |x| ^ oo) 
in which the peak amplitude becomes a two-valued function of the pulse duration. These 
solutions were proven to be stable as solutions to the CQNLS (@, [|^], |^], []T|). 

Equation (pTl| ) describes an idealized fiber; therefore, it is natural to consider the per¬ 
turbed CQNLS (PCQNLS) 

iAt = (1 -|- iea)Axx + iebA -|- (1 -|- iedi)\A\‘^A -|- (a -|- ied 2 )\A\^A, (1-2) 


where 0 < e <C 1 and the other parameters are real and of 0(1) (l^lj). The parameter a 
describes spectral filtering, b describes the linear gain or loss due to the fiber, and di and d 2 


describe the nonlinear gain or loss due to the fiber. Note that (1.2) is a well-defined PDE 
for e > 0 only if o > 0. 

Solitary wave solutions to (^) are found by setting 


A{x, t) = A{x)e 


—iujt 


(1.3) 
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and then finding heteroclinic and homoclinic solutions for the ODE 


(1 + iea)A!' + {—u> + ieb)A + (1 + iedi)\A\‘^A + (a + ied 2 )\A\‘^A = 0, 


(1.4) 


where ' = djdx. Equation (1.4) has been extensively studied by many authors (Q 
|]^|, |]I^, @1, |@, [^]). These papers have been concerned with finding various types 
of solutions, including fronts (kinks), bright solitary waves, and dark solitary waves. The 
methods employed have been both geometric ([p[|, || 8 |, p^, pTH, p^) and analytic (|p^, 

11, HD- 


Bright solitary waves exist when there are solutions to (^]^) which are homoclinic to 
|A| =0. When e = 0 and to > 0 the wave is given by the expression 


A^(x) = 


ioj 


1 + — (3 cosh(2Y^ x) 


13 = — -auj. 


(1.5) 


Since it is being assumed that a < 0, a restriction on (3 is that 0 < /3 < 1. 
expression for the wave exists even for e 


An analytic 

> 0 (ll^, @, @); however, it will not be given 
here. Eor the purposes of this paper it is enough to know that the wave exists for all e > 0. 

It was previously stated that the bright solitary wave is a stable solution to Q- 
However, recent numerical work by Soto-Crespo et al |^] suggests that this wave becomes 
an unstable solution to (0) for e nonzero. The numerics suggest that this instability arises 
from the presence of a real eigenvalue for the linearized problem moving out of the origin and 
into the right-half of the complex plane. The primary purpose of this paper is to determine 
if this is actually the case, and to determine possible stability/instability mechanisms. 

When discussing the stability of the bright solitary wave, one must locate the spectrum 
of the operator L found by linearizing ( |1.2D about the wave. The essential spectrum is easy 
to determine (Henry |lj]). When e = 0, it resides on the imaginary axis with |ImA| > cn, 
while for e > 0 it can be shown to be located in the left-half of the complex plane if a > 0 
and 6 < 0 (equation Q)- 

The location of the point spectrum is more problematic. It is known that when e = 0, 
zero is an eigenvalue of multiplicity four, and there are no other point eigenvalues (Weinstein 
1^1, iQ). For e / 0, two of these eigenvalues will remain at the origin, due to the spatial 
and rotational invariance of the PCQNLS, while the other two will generically move and be 
of 0(e). If either eigenvalue moves into the right-half plane, then the wave will be unstable. 
Unfortunately, one cannot conclude that if both eigenvalues move into the left-half plane, 
then the wave is stable. The reason is that it may be possible for eigenvalues to move out of 
the essential spectrum and into the right-half plane for e 7 ^ 0. This topic will be the focus 
of a future paper. 

In this paper a determination is made as to the location of the 0(e) eigenvalues for 
0 < e <C 1. In order to accomplish this task, it is necessary to perform detailed asymptotics 
for the Evans function, E(A), at A = 0. The Evans function is an analytic function whose 
zeros correspond to eigenvalues, with the order of the zero being the order of the eigenvalue 
(Q, [^1). Since the null-space of L is at least two-dimensional, E{0) = E'{0) = 0. Thus, 


when expanded about A = 0, the Evans function satisfies 

E{X) = E"(0)^ + E"'(0)^ + eW(0)^ + 0(A% 
with all the derivatives being real-valued. It turns out to be the case that 


( 1 . 6 ) 


E"(0) = Bie^ + 0(e^) 
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(equation (H), while if E''{Q) = 0, then 

E"'{{)) = B2e + 0{e^) 

(equation (|^) ). Furthermore, due to a result of Weinstein (@, M), = 0(1) 

(Corollary |2.3| ). Thus, if one can determine the signs of Hi, B 2 , and -£^^(0), then the 0(e) 
eigenvalues can be approximately located. 

Equation (Q) defines a four-dimensional ODE phase space. The quantity E”{0) is 
related to the manner in which the stable and unstable manifolds of = 0 intersect in this 
phase space, and a calculation of this quantity is similiar to the calculation which leads to 
the orientation index (Alexander and Jones Q, Q). The calculation of E'"{0) is a different 
matter, however. Using the ideas presented in Kapitula ||^, it is shown that E"'{Q) has 
a relationship with the projection of a certain function onto the null-space of the linear 
operator. In other words, 

Pf = BsE'"{0)A^, 

where P represents the projection onto the null-space, / is a particular function which 
measures the manner in which the stable and unstable manifolds intersect, is a certain 
basis function of the null-space of L, and H 3 > 0 is a constant of proportionality. The ideas 
leading to the calculation of £''"( 0 ) are generalized in an upcoming paper, as they are of 
interest in their own right. 

Eor the statement of the main theorems, set 


Ai = 


A*(x) dx, 


where A(x) is defined in (^), and let 


A24 = ^, a* = - —(<.A24 - j). 


Since the wave A{x) depends on w, so do the above constants. Asymptotic expansions for 
these constants are given in Appendix A. Note that 


dui-^d2 

duj-^2A 


> 0 


for 0 < /J < 1 (Proposition |5(^ . 


Theorem 1.1 Let 0 < e <C 1. Let 0 < /J < 1. Suppose that di = where d\ is given in 


Remark 1.5. Further suppose that a > 0, and that 


Set 


d 2 — -oia < 0 . 

O 


dujKd 2 / , 1 \ ^ n 

b = - {d 2 - -aa) < 0. 

Oujl'^24 o 


If 0 > b* > b, then there is one stable real eigenvalue and one real unstable eigenvalue, both 
of which are 0(e). However, if 0 > b > b*, then there are two stable real eigenvalues and 
zero unstable eigenvalues of 0{e). Furthermore, except for the double eigenvalue at zero, 
there are no other eigenvalues which are of 0{e). 
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Remark 1.2 If either a < 0 or b > 0, then the wave is unstable due to the presence of 
essential spectrum in the right-half plane. 


Remark 1.3 If a = d 2 = 0, then the wave is unstable. Thus, Ginzburg-Landau perturba¬ 
tions of the cubic NLS will only support unstable solitary waves. 


Remark 1.4 Since a > 0 and a < 0, the constant b* will be negative if and only if d 2 < 0. 
Thus, one can say that d 2 <t) is a minimal stability condition. 


Remark 1.5 The wave exists if di = d\, where 


dl = ^a- A 246 - Arf 2 ((i 2 - ^aa) + 0{e) 


(Corollary 5.1). 


OJ 


The constant d\ given in the above remark depends on uj, i.e., d\ = d\{a,h,d 2 ,w}). Let 
be a fixed parameter value, so that for uj ^ co* the wave does not exist without varying 


d\. As a consequence of Corollary 5.10 and the work of Kapitula and Maier-Paape 
has the following theorem concerning the existence of A^-pulse solutions to the PCQNLS. 


Theorem 1.6 Let the assumptions of Theorem \1.1\ be satisfied with d^ = d\{a,b,d 2 , 0 J*). 
For each N >2 there exists a bi-infinite sequence with 

lim = LO *, 

\k\^oo 


such that when to = co^ there is an N-pulse solution to (1^). If b < b*, then the N-pulse 
is unstable, and there exist at least N unstable eigenvalues. 


Remark 1.7 The interested reader should consult for a more complete description of 
the dynamics associated with (l.f)- 


The paper is organized in the following manner. In Section 2 the Evans function is 
constructed and its asymptotic behavior is determined. Sections 3 and 4 are devoted to 
deriving expressions for the various derivatives of the Evans function at A = 0. In Section 5 
the calculations are performed. Section 6 completes the argument leading to the two main 
theorems of this paper. 


2 Construction of the Evans function 

After the transformation A —> Ae~'^^^, the PCQNLS can be written in travelling wave 
coordinates (z = x — ct) as 

iA); = (1 + iea)Azz + icA^ + (—to + ieb^A + (1 + + (o; + ied 2 )|A|^A, (2.1) 
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where A is a complex-valued function of the variables {z,t) G R x R’*'. Upon setting 
A = Ai + iA 2 , and denoting A = (Ai, A 2 ), e becomes the system 

JA.t = (/2+£o>^)A^^-|-cJA2-t-(—CA;/2+£^«^)A.-|-(/2+£di'^)|ApA-|-(Q;/2+6d2</)|A|^A, (2.2) 


where I 2 is the 2x2 identity matrix and J is the skew-symmetric matrix 



The above system can be rewritten as 

J Aj = BA.ZZ T cJA.Z F(^A., io, , (2-3) 


where 

and 


Let A represent the bright solitary wave solution to ( |2.3D which is known to exist when 
c = 0 ([p3], @). When e = 0, A = {Rq,0)'^, where 


B — I 2 A cdJ 


F[A, uj, e) — (—(u /2 + €bJ)A {I 2 A ediT) | Ap A -|- {o'l2 + ed2'^) | A. 


2/ . 4a; 16 

^O(^) 1 j /-t O 1 (n \ > X Q 

1 + yl — (3 cos\i(2^u} z) S 

(|2^). Linearizing about the wave, the eigenvalue equation is given by 

RA" + DFA(A,o;,e)A = AJA, ' = 

dz 

i.e., —JLA = AA, where 


(2.4) 


(2.5) 


-JL = -J{Bd‘^ + DFAiA,uj,e)). 


( 2 . 6 ) 


A routine calculation shows that the essential spectrum for the operator —JL, hereafter 
referred to as ae{—JL), is given by 

ae{—JL) = {A G C : A = —e{arf‘ — 6) ± {if' -I- uj)i, rj G R} (2-7) 

(Henry |Q). Thus, for e > 0 the operator — JL is sectorial if a > 0, and the essential 
spectrum is in the left-half of the complex plane if 5 < 0. This observation leads to the 
following assumption, which is minimal if the solitary wave is to be stable. 


Assumption 2.1 The parameters a and b are sueh that a > 0 and b < 0. 


After setting Y = (A, A'), the eigenvalue equation (2.5) can be rewritten as the hrst- 
order system 

Y' = M{X,z)Y, (2.8) 


where M is the 4x4 block matrix 

M{X,z) = 


0 I 2 

B-^{XJ - DFA{A,uj,e)) 0 
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For A G = C\(Te(L) there exist complex analytic functions Y|(A, z) and Y“(A, z), i = 1, 2, 
which are solutions to and which satisfy 


1. jh^o |Y|(A,^)|=0, Y^(A,z) = (Yf AY|)(A,z)/O 

2. lim |Y“(A,z)|=0, Y“(A,z) = (Y]^ AY^)(A,z)/O 

2—> —CO 

(Q). The Evans function is given by 


E(A) = (Y]^ A Y 2 “ A Yi^ A Y|)(A, z), 


(2.9) 


and by Abel’s formula is independent of z. The Evans function is such that for A G fl it 
is zero if and only if A is an eigenvalue, with the order of the zero being the order of the 
eigenvalue (|jl|). Due to the invariances of the PCQNLS, two solutions to when A = 0 
are A = A' and A = JA. As such, one can set 

1. Yf(0,z) = Y3^(0,z) = U' 

2. Y|(0,z)=Y2“(0,z) = Uj, 

where 

U = (A,A')'^, Uj = (JA, JA')^. 

The below lemma describes the asymptotic behavior of the Evans function. 


( 2 . 10 ) 

( 2 . 11 ) 


Lemma 2.2 

as A —> oo. 


With Y“(0, z) and Y|( 0 , 2 ;) as described in ( 2.1C ), i/A G R, then E{\) < 0 


Proof: It can be assumed without loss of generality that e = 0, as if the result is true for 
e = 0, it will then be true for 0 < e <C 1. Assume that A G R^. 

Let Y = (P, Q)^ in (^). Upon setting s = \/Az and Q = \/AQ and letting A —> oo, 
equation (| 2 . 8 D becomes the autonomous system 


P 

Q 


0 h](P\ 

J 0 J Q J ’ 


where ' = d/ds. The eigenvalues of the above matrix are given by 7 (±l±z), where 7 = \/2/2, 
so that there exists a two-dimensional unstable subspace and two-dimensional stable sub¬ 
space, with the two-dimensional unstable subspace being Span{( 7 , — 7 ,1,0)^, ( 7 , 7 , 0 ,1)^} 
and the two-dimensional stable subspace being Span{( 7 , — 7 , —1,0)^, (— 7 , — 7 ,0,1)^}. 

Let Oj A ej = eij. In A^(R'^) the unstable subspace is represented by the vector 

Y“(-hoo) = ( 7 ,- 7 ,1,0)^ A ( 7 , 7 ,0,1)"^ 

= ei2 — 7613 -I- 7614 — 7623 — 7624 -I- 634, 


while the stable subspace is represented by the vector 


Y^+oo) = ( 7 ,- 7 ,-1,0)^ A (- 7 ,- 7 ,0,1)^ 

= -ei2 — 7^13 + 7^14 — 7^23 — 7^24 — 634. 


Note that 

Y“(+oo) A Y"(+oo) = -2. 


( 2 . 12 ) 


6 





When A = 0, for each fixed z both the unstable and stable subspaces are spanned by 
the vectors (i?Q( 2 ;), 0, i?Q(0), 0)^ and (0, 0,-Ro(-z))^. Set 


Y“(0) = lim e“2^^Y“(0,z) 

2—OO 

Y^(0) = lim e2^^^Y^(0,z). 

2—>00 


(2.13) 


Using the representation for Rq it can then be seen that 

Y“(0) = lim A{0,Ro{z),0,R'o{z)f 

2—^ —CO 

= 0,w,0)'^ A (0,1,0, \/a;)'^ 

= ^(^ei2 + a;ei4 - a;e23 + w^/^e34), 

where 

8 ,, 16 

A similiar calculation shows that 

Y^(0) = lim e2V^"(i2(,(z),0,i?(,'(0),0)^ A(0,i2o(^),0,i2()(^))'^ 

2—>CO 

= //(-Vwei2 + a;ei4 - a;e23 - w^/^e34). 

Note that 

Y“(0) A Y"(0) = -4a;V^- 

A more complete discussion of the following argument can be found in Alexander and 
Jones Q. An orientation of is given by a nonzero element r/ of A^(R^). Two ordered 
bases of R^, {wi, ..., W4} and {ui,..., U4}, have the same orientation if wi A • • • A W4 is a 
positive multiple of ui A • • • A U 4 (Q). Any two such bases are related by a matrix having 
positive determinant. 

The functions Y“(A, 2 ) and Yf(A, 2 ;) can be used to get a basis for R^ for each A and 
z. In particular, in a manner similiar to (|2.13| ) proper scalings of 

lim Y“(A, 2 ;), limY|(A,z) 

z^—00 2—>00 


determine a basis for each A. Since 


[Y“(0) A Y"(0)][Y“(+oo) A Y"(+oo)] > 0, 

the matrix taking the basis at A = 0 to that at A = +00 has positive determinant, so that 
the two bases have the same orientation. The sign of E{X) for large positive A is then 
determined by equation ( |2.12| ), from which the conclusion of the lemma follows. ■ 


Corollary 2.3 When e = 0, the Evans function satisfies E{0) = E'{0) = E"{0) = E'"{0) = 
0, with < 0. Furthermore, E{X) < 0 for A > 0. 


Proof: The fact that the first three derivatives of the Evans function at A = 0 are zero, 
with the fourth derivative being nonzero, is a direct consequence of the work of Weinstein 
( p^ , 1^^). Furthermore, it is known that when e = 0 the bright solitary wave is stable, so 
that there exist no positive eigenvalues (|^, |jl^); hence, the Evans function is nonzero for 


A > 0. The fact that the fourth derivative is negative then follows from Lemma 2.2 
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3 Calculation of derivatives 


For e > 0 it will be generically true that -E(O) = E'{0) = 0 with E''{0) / 0. Since 
E^^^{0) < 0, by calculating E''{0) one will be able to determine the location of the zeros of 
E{X) which are 0(e), and hence the location of the small eigenvalues. When E”{0) = 0, an 
eigenvalue will be passing through the origin. A determination of E'”{0) will enable one to 
decide whether the eigenvalue is passing into the right-half or left-half of the complex plane. 
This section is devoted to determining these quantities, and relating them to properties of 
the wave. 

Time independent solutions to (|2.3|) satisfy the ODE 

fiA" + cJA' + F(A,w,e) =0, ' = (3-1) 

which can be written as the first-order system 

U'= G(U,c,a;,e), (3.2) 

where U = (Ui,U 2 ) G and 

G(U,c,(u,e) = ) ' 


The bright solitary wave corresponds to a solution homoclinic to U = 0, and is realized as 
the nontrivial intersection of the two-dimensional unstable manifold, W'^{z, c, uj, e), with the 
two-dimensional stable manifold, W^{z,c,uj,e). Due to the rotational symmetry associated 
with the PCQNLS, there exists no distinguished trajectory in W^lz, c, uj, e) n W^(z, c, uj, e). 
However, this rotational symmetry allows one to choose a trajectory so that ^ 2 ( 0 ) = 0, 
which uniquely defines a trajectory in the two-dimensional manifold. Set U = (A, A'), so 
that U C W^(z, c, UJ, e) n W^{z, c, uj, e) is a distinguished solution. 

Before continuing, the following proposition is needed. It follows immediately upon 
examination of (2.3). 


Proposition 3.1 The Frechet derivative of the nonlinearity F satisfies 

DF^(A,a;,e) = - A. 


(3.3) 


Since G depends smoothly on the parameters, so do the manifolds. The bright solitary 
wave is manifested as the nontrivial intersection of W'^{z,c,uj,e) and VF®(z, c,a;, e). Dif¬ 
ferentiating ( p.2|) with respect to the parameters c and uj and evaluating over the wave U 
yields the systems 


{dcW^y = DGu{A, 0, UJ, e) dcW'" + (0, -B-^JA'f (3.4) 


and 

{d^W^y = DGu{A, 0, w, e) 5^ W" + (0, B-^Af. (3.5) 

In these equations r G {u, s}, the result of Proposition ^T| is implicitly used, and 


DG\j{A,0,uj,e) 


0 I 2 

—B~^DFx{A,uj,e) 0 





Note that a consequence of these equations is that dciW'^ — W^) and — W^) are 

solutions to the linear system 

5U'= L>Gu(A,0,a;,e)(5U. (3.6) 

If one sets Uj = (JA, JA'), then the following proposition is realized. 

Proposition 3.2 Four solutions to are given bytl', Uj, and d^{W^- 

W^); furthermore, if 

D 2 = [dc{W^ -W^) Ad^iW'^ -W^) AJj' Aijj]{z,0,u;,e) 
is nonzero, then the solutions are linearly independent. 

Proof: It has already been seen that these four functions are solutions to the linear system. 
When D 2 0, the linear independence of the solutions follows from the fact that D 2 is the 
Wronskian. ■ 


Remark 3.3 By Abel’s formula, D 2 is independent of z. 

Set 

= dc{W^ - W^), 6 V 2 = - W^), 6 V 3 = U', 6 V 4 = Uj, (3.7) 

so that 

L>2 = (<5Ui AhUa AdUs A(IU4)(z,0,w,e). 

Assuming that D 2 7 ^ 0, the functions 6 Ui and 6 U 2 grow exponentially fast in the supremum 
norm as \z\ —> 00 , while the functions 6 U 5 and 6 U 4 decay exponentially fast. 

Let H : R ^ be a uniformly bounded measurable function. Suppose that the 
solution to 

hU'= L»Gu(A,0,u;,e)<5U + H (3.8) 

is desired, and further suppose that one wishes the solution to be bounded for either z —> 
—00 or z —> 00 . Denoting the solution by with | 5 U^( 2 ;)| < M < 00 as 2 ; —> ± 00 , by 

following the discussion in Kapitula it can be seen that 


1 


= —(c±(H)hUi + cf (H)hU 2 + C 3 (H)hU 3 + C4(H)hU4), (3.9) 

F>2 


where 

cfm = 


|HhU 2 5U3hU4|(s)ds, c^(H) = 


±00 


/ ±00 


IhUiHhUg 6V4\is)ds 


C 3 (H) = / |5Ui<5U2H(5U4|(s)ds, C 4 (H) = / |(5Ui JUa (5U3 H|(s) ds. 

Jo Jo 

(3.10) 

The following lemma can now be proved. 


Lemma 3.4 Set 




Hi 

= {0,-B-^JA' 

)^, 

Then 


POO 




D2 = 

j —00 

dUa (5U3 dU4|(s: 

) ds 

while 


POO 




0 = 

/ IH 2 

J —00 

hUa dUs hU4|(s) 

ds 


A ^T 
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Proof: Using equations (^) and along with equation (|3.9|) , one can see that 
dcW'^ = ^(cr(Hi)(^Ui + C2(Hi),5U2+C3(Hi)5U3 + C4(Hi)5U4) 

D2 

dcW^ = —(c+(Hi)(^Ui + c+(Hi),5U2 + C3(Hi) 5U3 + C4(Hi)JU4), 

L>2 

and 

= ^(cr(H2)(^Ui+C2(H2)<5U2+C3(H2)<5U3 + C4(H2)5U4) 

d^W'^ = —(c+(H2)(^Ui + C+(H2)<5U2 + C3(H2)<5U3 + C4(H2)5U4). 

L>2 

Subtracting and using the definition of (5Ui and (5U2 then yields 

5Ui = ^[(cr(Hi)-c+(Hi))(lUi + (c2(Hi)-c+(Hi)),5U2] 

D2 

5 U2 = —[(cr(H 2 )-c+(H 2 ))(lUi + (c2(H2 )-c+(H2 )).5U2]. 

Using the definitions of the c^’s and the fact that (5Ui are linearly independent functions 
yields the final result. ■ 


Define 

e* =-JUi A5U3 A(5U4, 62 =-(IU 2 A(5U3 A<5U4. (3.11) 

The functions e* £ A^(R^) for i = 1,2; furthermore, since A = 0 is an isolated eigenvalue, 
both of these functions satisfy an estimate of the type 

|e-(z)| < i = l,2 (3.12) 


for some positive constants C and jjL (Kapitula j^). For a given bounded continuous 
function F : R ^ R^, define 

/ CO 

(HAe-)(s)ds, i = l,2, 

-00 


where H = (0, F)'^. 
following manner. 


With the above discussion in mind, one can rewrite Lemma 3.4 in the 


Corollary 3.5 The constant D 2 is given by 

D 2 =< el,B-^A >=< e*2,B-^JA’ > . 

Furthermore, 

0 =< ei,R“^JA' >=< e*2,B-^A > . 

It is now possible to relate the Evans function to the structural stability of the wave. 
The proof of the first part of the below lemma is an alternate to that found in Alexander 
and Jones |Q, and may be of interest in its own right. 

Lemma 3.6 The Evans function satisfies 

E"( 0 ) = 2 D 2 , 

where D 2 is defined in Proposition \3 .4 Alternatively, 

E"{0) = 2 < el,B-^A >= 2 < 62 , B-^JA' > . 
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Proof: In this proof, the dependence of functions on the variable z will be supressed. Upon 
differentiating E{X) and evaluating at A = 0, one sees that 

U"(0) = 2(aA(Yl‘ - Y^) A 5 a(Y^ - Y|) A Y? A Y|)(0). 

In the above calculation the fact that Y“(0) = Y|(0) was implicitly used. Since Yf (0) = U' 
and YKO) = Uj, all that is left to do is show the equivalence with the first two entries 
making up E"{0). 

Differentiating ( p.8D with respect to A and evaluating at A = 0 one sees that 
(5 aY[)' = M(0, z) dxYl - (0, 


and 

(5aY^)' = M(0, z) 5 aY2 ^ - (0, 


where r G {u,s}. In the above equation, the fact that = —I 2 is implicitly used. Since 
M(0) = iAGu(A, 0, a;, e), by following the proof of Lemma ^ and using the definitions of 
Hj presented therein it can be shown that 


5A(Y3^-Yf) 

5a(Y2“-Y|) 


-^[(cr(Hi) 

D2 


-c+(Hi))(IUi + (c2(Hi) 
-C+(H2))(^Ui + (c2(H2) 


-C+(Hi))5U2] 

-C+(H2))5U2]. 


By the result of Lemma [}.4| it is then seen that actually 


5a(Y“ - Yf) = -,5Ui, 5 a(Y2“ - Y|) = -,5U2. 


Upon substituting the above into the expression for U"(0) the first part of the lemma is 
proved. 

The second part of the lemma follows immediately from Corollary ■ 


Remark 3.7 Note that a consequence of the above argument is that 

dxY{ = -5cIU^ 5a Y^ = -5^IU'’ 


for r G {u, s}. 


Remark 3.8 Using the expansion of the Evans function given in (1^) one can then write 


E{X) = <el,B-^A> X^ + E"'{0)^ + E^^\0)^ + O{X^) 


= < e^, B-^JA' > A2 + E"'(0)^ + E(^^(0)^ + O(X^). 


The above discussion is predicated on the assumption that B 2 / 0, which is equivalent 
to the manifolds and VU® intersecting transversely. There are instances in which this 
intersection will not be transverse, in which case E >2 = 0. In this circumstance E"(0) = 0 


(Lemma 3.6), so that an eigenvalue is passing through the origin. In order to determine the 
direction in which the eigenvalue is moving through the origin, it would be helpful to know 
E"'(0). 
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Suppose that D 2 = 0 due to the fact that 5 U 2 = 0, i.e., because , while 

e\ / 0. Note that this implies that —> 0 exponentially fast as \z\ 00 . Let (5U2 be 

any solution to ( |3.6|) such that 

Ds = { 6 \J 2 /\ el){z,0,u;,e) (3.13) 

is nonzero. The above discussion concerning concerning the construction of solutions to 
(|3.8|) can be recreated by substituting < 5 X 12 with ( 5 U 2 and D 2 with ZI 3 . Let tt : be 

the projection operator onto the first two components. The following lemma can now be 
proven. 


Lemma 3.9 Suppose that E"{0) = 0 with el 7 ^ 0. Then 

E"'{0) = 6 < el,B-^j7r{d^W^) > . 

Proof: In this proof the dependence of solutions on z will be supressed. Since E"(0) = 0 
and 63 7 ^ 0 implies that ( 5 U 2 = 0, by the proof of Lemma it is necessarily true that 
d\C ^2 ~ Y 2 )( 0 ) = 0. Using this fact, a tedious calculation then shows that 

E"'{0) = 3{dx{Y^ - Yl) A - Y|) A Y^* A Y|)(0). 


Define the projection matrix 


Q = 


0 0 

B-^J 0 


and set H 3 = Qd^jW^ = {0,B ^ j 7 r{dcoW'^))'^. Note that Q = M\{0,z). By Remark ^ 
Using this, differentiating 


9a Y| = 9a Y“ = -9^1UT 

twice with respect to A, and evaluating at A = 0 gives 

{dlY^^y = M(0, z)dlY^^ - 2Q9^1U^ 

where r £ {u, s}. Using the definition of H 3 , the solutions to this ODE are 

92Y|(0) = - A(c-(H3)5Ui + C 2 (H3)5U2 + C3(H3)<5U3 + C4(H3)(5U4) 

T>3 


dlYm = -^(c+(H3)5Ui + C+(H3)5U2 + C3(H3)<5U3 + C4(H3)<5U4), 

where the above functions c* are such that in their definitions (5U2 has been replace with 
(5U2. Upon subtracting one gets that 

9|(Y2“ - Y|)(0) = -^[(cr(H3) - C+(H 3 ))<^Ui + (c2-(H3) - C+(H3))5U2]. (3.14) 

By the previous lemma it is known that 

9A(Yj^-Yf)(0) = -,5Ui. 

After substituting the above expressions into that for E'”{0) and using the definition of 
one gets that 

E'"(0) = 6(c2-(H3)-c+(H3)). 

Evaluating this expression gives the final step in the proof. ■ 

Remark 3.10 Using the expansion of the Evans function given in if E"{0) = 0, then 

E(A) =< el,B-^j7r(d^W^) > + e(^^(0)^ + O(A^). 
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4 Alternative expressions for the derivatives 


Now that expressions are known for the various derivatives of the Evans function, it is 
desirable to reduce them to computable quantities. This section is devoted to that task. 

It will be convenient to write everything in polar coordinates, i.e.. 


A = (r cos 0, r sin 9). 


In polar coordinates, U = T{r, 9, s, (p), where 


T{r,9,sA) 


r cos 9 
r sin 9 

rs cos 9 — rcj) sin 9 ’ 

rs sin 9 + rcj) cos 9 


with s = r'/r and (/> = A routine calculation shows that 


DT{r,9,s,^) 


cos 9 —r sin 9 

sin 9 r cos 9 

s cos 0 — (/) sin 0 —r(s sin 0 + (/> cos 0 ) 
s sin 0 + (/) cos 0 r(scos 0 — (/)sin 0 ) 


0 0 

0 0 

r cos 9 —r sin 9 ’ 

r sin 9 r cos 9 


(4.1) 


with 

\DTir,9,sA)\=r^ 

so that the transformation is nonsingular except at the origin. 

In polar coordinates, let the manifolds be denoted by and W^. In these coordinates 
it is a routine calculation to show that 


,5Ui = DTdc{w;^-wp, dUs = DTd,W}^ 

6\J2 = DTdUW^-Wp, 6Ui = DTdeW^. 

When s = 0 the manifolds can be parameterized as 

W; = {r^{9,(5),9,0,cp^{9,l3)f, w; = {r^{9, (3),9,0,cP^{9, P)f, (4.3) 

where (5 = (c, cu, e). Now, let the underlying wave be denoted by 

A(z) = (/^(z) cos ©(z), i?( 2 ;) sin ©(z)). 

Due to the fact that the wave is even it can be assumed that 

i?'(0) = ©'(0) = 0, (4.4) 

while the rotational symmetry of the PCQNLS allows one to set 

©(0) = 0. (4.5) 

Under these assumptions, when z = 0, i.e., when S = R'/R = 0, 

dc{w; - w;) = dc{{r^ - r^), 0,0, (0“ - cl>^)Y 

5^(Wp“ - w;) = 5^((r“ - r^), 0,0, (4.6) 

a,IT“ = (0,0,5'(0),<I>'(0))^, deW;^ = (0,1,0,0)^. 

Combining the above with (|4.2|) allows one to prove the following lemma. 
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Lemma 4.1 The Evans function satisfies 

E"(0) = -8ii2(0)ii"(0)5cr“(0,0,u;,e)5a;</'“(0,0,u;,e). 

Proof: First note that 


L >2 = {\DT\ dc{W^ - wp A dpw^ - wp A A 5,VF“)(0,0, o;, e). 
By (|4.6| ) and the calculation for \DT\, it can be seen that 


L>2 = -R^{0)S'i0) 


dc{r'^ — r*) dPr'^ — r^) 


The steady-state equations in polar coordinates are given by equation (|5.3|). As a con¬ 


sequence of Proposition 5.1 


r“(0, c, O', e) = r^(0, —c, cu, e), c, w, e) = ——c, w, e). 

Thus, it can be concluded that 

dPr'" - r^)(0,0,a;,e) = dfifT- - (/)®)(0, 0, w, e) = 0, 

while 


<9c(r“ - r^)(0,0,a;,e) = 2(9cr“(0,0, w, e), - (/>'')(0,0, w, e) = 2cl^(/)“(0,0, w, e). 

Since 5''(0) = ii"(0)/ii(0), this then yields that 

D2 = -4i?2(0)i?"(0)acr“(0,0,a;,e)a^<^“(0,0,a;,e). 

The statement of Lemma |3.6| then gives the result. ■ 

Now that a computable expression for E" {0) is known, it would be beneficial to have an 
expression for E'''{fi). Observation of Lemma |3.9| yields that one must first better understand 

e* = -SUi A SUs A SU 4 e A3(R^). 

By p.2| ) the above quantity can be rewritten as 

et = -(DTdc(Wp--Wp)A(DTd,WpA(DTdeWp 

= -DT(^pdfiWfi-Wp Ad^W^ AdffWfi), ^ ^ 

where is the 4x4 matrix induced by DT which maps A^(R^) to itself. The matrix 

is formed by taking all the 3x3 minors of DT, and is given by 

(s cos 0 -I- (/) sin 0) — {s sin 6 — 4> cos 9) 

s sin 9 — (f> cos 9 — (s cos 9 + 4) sin 9) 

cos 9 sin 9 

—r sin 9 r cos 9 

(Q, 1^). Thus, in order to finish the calculation of e\, all that is left to determine is 

{el)p = 9c(VF; - Wp A d,W;f A deWp 


jjTi^) = 


cos 9 sin 9 
— sin 9 cos 9 
0 0 

0 0 
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Set 


^i = d,w; = {R',Q',s\^'Y 

^2=dcW^, ^t = dcW; (4.8) 

6 = 9eW'“ = (0,l,0,0). 

Let the vectors ej, i = 1,... 4, be the unit vectors in R'^, and define 


^ijk — ^ ®j • 

The collection of vectors { 0123 , 6124 , 6134 , 6234 } form a basis for A^(R^), so that (e{)p can 
be written in terms of these vectors. Now define 


p±^ i^l)^ i^t)^ 

( 6 ), (^ 2 ^), 


(4.9) 


and set Pij = P^j — Pj^. Using (|4.8|) , a routine calculation then shows that 

(ei)p = ^3^123 + -P146124 — -P346234- 
A consequence of the above discussion is the following lemma. 


Lemma 4.2 Let M > 0 be given, and suppose that Pij = 0(e) for \z\ < M. Then 

* _ i —R^{Pi3^123 + (-P 14 + 5'.P34)6i24 — RP3i^23i) + 0(e^), \z\ < M 

“ 1 0(e-^l^l)e, |z| > M, 


where /U > 0 . 


Proof: Let 7/1 > 0 be such that R^(z) = 0{e i.e., r/i = 2y/uj + 0{e). Since = 0(e) 

implies that 0 = 0(e)z, one can easily see that for \z\ < M 


Dt(^) = i ?2 


10-50 
0 10-5 
0 0 1 0 
0 0 0 R 


+ 


while ||Dr(^)|| = 0(e '^il^l) for \z\ > M. Therefore, given the assumption on the functions 
Pij, for \z\ < M, 

et = Zir(3)(e|)p_ 

= A36123 + (A4 + 5 P 34 )ei 24 — P-P346234 + O ( e ^). 

Since Pij = 0(e) for \z\ < M, it is necessarily true that for \z\ > M, Pij = 0(e'^^^^^)e, 
which then implies that (e|)p = 0(e’^^l^l)e. Thus, for \z\ > M one sees that 

|e}| < ||Pr(3)|||(e})p| 

= 0(e-’?il^l)0(e^2|d)e 

= 0 (e(’? 2 -m)ld)e. 

Setting n = rp — ri 2 , the fact that e} approaches zero exponentially fast (equation (|3.12|) ) 
guarantees that p > 0. ■ 
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Since B = I 2 + eaJ, a simple calculation shows that B ^ J = J + 0(e). Thus, using the 
fact that 0 = 0 (e) for \z\ < M, it is not difficult to see that when ( 5 U 2 = 0 , 




{0,d^Ro)^+ 0{e), \z\<M 


|z| > M, 


(4.10) 


where 773 > 0 . 


Lemma 4.3 Let M > 0 be given, and suppose that Pij = 0(e) for \z\ < M. When 
E”{0) = 0, the third derivative of the Evans function satisfies 

/ OO 

Rl{s)dMs)Pi3{s) ds + 0{e-^^^)e + 0{e^), 

-OO 

where 7/4 > 0 . 

Proof: The integrand associated with £''^^(0) is given by H 3 A e^, where 

H3 = iO,B-^Jn{d^W^)f. 


Using ( 4 . 1 c) and Lemma 4.2, one then sees that for \z\ < M 

H 3 A 64 = —^4123 + O(e^) 

= Rld^RoPi3 + 0{e^), 

while for \z\ > M 

Us A el = 0(e-^4^l^l)e. 

In the above calculation, the fact that 64423 = ©4 A 6423 = — 1 is used. 
Remark 4.4 A similiar calculation leads to the conclusion that 

roo 

£"(0) 2 / Ro{s){Pu{s) + So(s)P34(s)) ds. 


5 Asymptotics 


Now that an expression for £"(0) has been derived, in order to determine the location of 
the eigenvalues near zero the expressions dcr^ and must be calculated. In addition, 

in order to calculate £'^^( 0 ), one must determine Pij and show that the quantities are 0 (e) 
for \z\ < M. 

Set 

A(z) = {r(z) cos 9{z),r{z) sin 6(z)). 

Let the known underlying solitary wave be denoted by (i?, 0, S, <h). Note that 


S = R'/R = -^InR. 
dz 


(5.1) 


Recall the analytic expression for the wave given in (2.4) when e = 0, i.e., 

4a; 


Rn(z) = 


Qo{z) = 0 . 


_ g ^ _ 16^^ 

1 + \/l — j3 cosh(2\/a; z) ’ 3 


(5.2) 
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This wave will henceforth be denoted by (i? 0 ) 0, ^o, 0)^. By defining 


A = l + e2a^ 


the steady-state ODE is 

r' = rs 
6' = (j) 

As' = —As^ + A0^ — c{eas — (f>) — {—cv -|- e^o6) (5-3) 

— (1 + e^adi)r^ — {a + e^ad 2 )r^ 

A(j)' = —2As4> — c{s + eacj)) — e{{b + aoj) + {di — a)r^ + {d 2 — aa)r^]. 

It should be noted that the equation for 6 is superfluous, and hence is usually ignored; 
however, it is included here for completeness. After dropping the 0{e^) terms the variational 
equations are given by 

6r' = S6r + R6s 

66' = 6(1) 

6s' = —2R{1 + 2aR‘^)6r — 2S6s + 2^6(j) + 6u) — {eaS — ^)6c 

—2eo[(6 -|- auj) + (di — a)i?^ + (^2 — aa)R‘^]6€ 

6(j)' = —2eR[{di — a) -|- 2(^2 — aa)R^]6r — 2<h(5s — 2S6(j) (5-4) 

—ea6uj — (ea<I> + S)6c — [(6 + aoj) + (di — a)R? + (^2 — a(y)R‘^\6e 

6u)' = 0 

(5c' = 0 

6e' = 0 . 


An observation yields the following proposition. 


Proposition 5.1 Equation (5A) is invariant under 

{z, c, w, r, (9, s, (j)) {-z, -c, cu, r, 6, -s, -cp). 

An expression for will first be determined. Let 

(j)^{z) = 6(j){deWp{z,c,uj)). 


(5.5) 


Since 

6c{deWp{z,c,u))) = 6uj{deWp{z,c,u))) = 0 , 

by using (|5.4D it can be seen that when e = 0 

4>'^ = -2SQct)e - [{b + auj) + {di - a)Rl + (^2 - aa)RQ). (5.6) 

By definition (/>£ is uniformly bounded as 2 —> —oo, so that upon using (P) the solution to 
(|5.6|) can be written as 


Rl{z)(t>e{z) 


-[{b -I- auj) [ 

J —OO 

+ ((^2 - 0 . 0 ) 


^^ 0 ( 5 ) ds “h {d\ — 

r Ri{s)ds]. 

j —00 



Rq{s) ds 


(5.7) 
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By definition the function cj)^ describes, up to 0(e), the location of the (^-component of 
Wp{z,c,ijj), so that 


'/’“(O) = e(/>,( 0 ) + 0{t^). 


(5.8) 


Thus, when performing calculations on (/>“(0), for e > 0 small enough it is sufficient to 


perform them on i?i>e( 0 ). Given (^) and the fact that an exact expression exists for i?g(z). 


this then implies that rather detailed information can be gathered regarding the variation 
of </>“(0) with respect to u) for e sufficiently small. Set 



(5.9) 


Lemma 5.2 The function is given by 


(/>“( 0 ) — e(/>e( 0 ) + O(e^), 


where 


2 i?g( 0 )(/)£( 0 ) — A 2 O — A 26 — h.^d\ — Ag(i 2 - 


Proof: Since Rq is an even function. 



for any positive integer m. Thus, when (^]^) is evaluated at z = 0, 

2i?g(0)(/>£(0) = {—u)h.2 + A 4 + aAg)a — A 26 — K^di — Ag(i 2 - 


The function Rq satisfies 


R'q - wRq + + aRl = 0 . 


Upon multiplying the above equation by Rq and integrating by parts one sees that 


A2 — —CUA2 + A4 + oAg, 


from which the conclusion of the lemma follows. ■ 

Remark 5.3 The expressions A 2 and A 4 are evaluated in Appendix A. 

It is known that the wave exists for all e > 0, with the perturbation being regular 
(p^, [p5| ). A necessary condition for the existence of the bright solitary wave is that <fe{z) 
remains uniformly bounded as 2 ; —> 00 . Since Rq{z) —> 0 as z —> 00 , this then yields the 
next lemma. 

Lemma 5.4 A necessary condition for the existence of the bright solitary wave is that 


A^fl — A26 — A 4 (ii — Ag(i2 — 0- 


Proof: Evaluating ( |5.7D at z = 00 and requiring that the right-hand side be zero at the 
limit yields 


0 — (—cuA 2 -|- A4 -|- Q;Ag)a — A26 — A 4 ^di — Ag(i2 
= A^ffl — A26 — A^di — Ag(i2- ■ 
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Remark 5.5 An examination of Lemmas 5., 
for the existence of the wave implies that 


and 5.4 shows that the necessary condition 


4(0) = 0. 


The expression present in the above lemma can clearly be solved for di in terms of 
the other parameters. Before doing so, however, it will be desirable to simplify the above 
expression. As the following proposition illustrates, there is a simple relationship between 
the above quantities. 


Proposition 5.6 The relations 

1. Ag 

2. A '2 


hold true. 





Proof: As mentioned in the proof of Lemma the function Rq satisfies the ODE 

Rq - coRo + Rq + aRl = 0. 

Multiplying by Rq and integrating by parts yields that 

—A 2 — ujA.2 T .^4 T oAg = 0, 
while multiplying by R!q and integrating yields 

1 1 

A2 — u}K 2 + —A4 + —oAg = 0. 

Upon subtracting the above two equations one sees that 


./ 1 . 1 . 

A 2 - -A 4 - -oAg = 0. 

The conclusion of the first part of the proposition is now clear. 

The proof for the second part follows in a similiar manner. Simply add the two equations 
to get the relation 

3 2 

U)K2 — ^^4 — -oAg = 0, 

from which one immediately gets the second part of the proposition. ■ 

Note that for fd = — IGaw/S the relation for Ag can be rewritten as 

Sw 3 

Ag - —- 4 ^ 4 )- 

With this observation, dehne 

Am = h, Aj, = ~{uAm - j). (5.10) 

Note that as a consequence of Proposition |5.6|, A^j = Ag/A 4 . 
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Corollary 5.7 In order for the wave to exist the parameter di must equal d\, where up to 
0 {e) ^ ^ 

dl = -a- A 246 - Arf 2 (d 2 - -^aa). 


Proof: By Lemma 5.4, in order for the wave to exist it must be true that 

, A2 A2 Ag 

“1 = - 7 — 02 - 

A 4 A 4 A 4 

The conclusion of the corollary follows after one uses the relationships described in Propo¬ 


sition 5.6 


Now that an expression for the function (/>“(0) is known (Lemma 5.2), it is possible to 
understand its behavior when the parameter w is varied. A consequence of Lemma is 
that it is sufficient to understand the manner in which 4>e{0) varies. Since </>e(0) = 0 when 
d = di, a simple application of the implicit function theorem yields that 

dujf'eiO) -h dd^(j)eiO)dujdl = 0. (5.11) 

The quantities dd^cfeiO) and d^^dl are accessible, so that the term cl^(/)e(0) can be calculated. 
Lemma 5.8 When di = d\, 


duj(t>e{^) = 


A 4 


2Rl{0) 

Proof: First, an examination of Lemma 


ddi4>e{0) = - 


dui^- 2 ib + d^Ad 2 {d 2 - -aa) 


shows that 
A 4 


2 i? 2 (o)' 


Upon differentiating the expression for d\ given in Corollary 5.7 and using (5.11), one then 
arrives at the conclusion of the lemma. ■ 

It is important to understand how 9tj(^e(0) varies with the parameters. Before making 
a definitive statement, the following proposition is needed. 


Proposition 5.9 Set 

When 0 < P < 1, 
so that 


P = — -au). 


5<^A24 < 0, d^kd^ > 0, 

duiAd2 ^ Q 


dwA24 

Proof: Using the definition of A 24 and Corollary |A.2| one sees that 

(A 4 — 4a;A2)(9a;A2 


9a;A24 = 


AI 
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Upon using the Taylor expansions given in Corollary [A.3| one sees that 

(2n + l)(2n + 3)’ 


A 4 — 4a;A 2 = — 16a;^'^^ ^ 


which is clearly negative. Since Corollary A.2 states that d^^K 2 > 0 for 0 < /9 < 1, it is now 
clear that < 0 . 

Since = /3/a;, 

du}-^d2 — ~^dij(^UjA.24.') ■ 


Using the Taylor series expansions given in Corollary [A.3| , after some tedious manipulations 
one can see that 


dUioA2i) = C^{an-bn)l3^, 


where 


and 


n=0 


1 


-1 


>0 


n=0 


H 9.- 


j '^0 2/ + 1 2(n - j) + 3 


’ 9.- 


^0 2j + 3 2(n - j) + 1 ■ 


The claim regarding will be proven as soon as it can be shown that a„ — < 0. 

Upon combining terms, 


— 4^ 


. n-2j 


where 


f{j, n) = {2j + l)(2(n - j) + l)(2j + 3)(2(n - j) + 3). 
By the integral test, 

rn 

CLn — bn < 4: / xg{x, n) dx, 

Jo 

where 


gix,n) = 


n — 2x 
fix,n)' 


Set y = X — n/2. Then 


g{y,n) = -2 


fiy,n) 


with 


1 


/(y> n) = A( 4?/2 - (n + l)^)(4y2 - (n + 3)^), 


so that g{y, n) is odd in y with yg{y, n) < 0. Therefore, 

rnl2 
l-nl2 

rn/2 

J—nl2 

< 0 , 


rn /'^/2 Yl 

/ xg{x,n)dx = / (y + -)9iy^'^) dy 

Jo J-n/2 2 

ygiy, n) dy 
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so that ttn — bn < 0 . 


Combining the above results yields the following corollary, which concerns the variation 
of with (jj. 

Corollary 5.10 Suppose that d = dl- Set 

duj^d2 ( 1 1 \ 

^ “ 5““’- 

For e > 0 sufficiently small, if b > b*, then > 0; otherwise, < 0. Further¬ 

more, for 0 < /3 < 1 

duj-^d2 


dnjh.' 


< 0 . 


Proof: Since A 4 /i 2 g( 0 ) > 0 and 9ijA24 < 0, the result follows immediately from Lemma 
5.8 and Proposition ■ 

Now that 9aj0“(O) is known, the quantities 5cr“(0) and Pij must be calculated. This 
can be accomplished simultaneously. As in (O), set 

6 = ff = d,w;f, ct = dcW;. 

Letting Px^xj denote 6xi A Sxj, as in ([4.91) set 

Pri=Prs{Clffit), Pr% = Prcti^l,Ct), P^^ = PsctiClffi^) ■ 

Note that the computation of E'”{0) requires that Ppg — P+ be known (Lemma ^) . Before 
continuing, a preliminary lemma is needed. 


Lemma 5.11 Set 


When e = 0, (ff = —1/2 


= mt 


Proof: It is easy to see from the variational equation (5.4) that when e = 0 

{f)fy = -2Soff - Sq. 

This equation is easily solved, and one then finds that 


(5.12) 


Poiz)(ftiz) = -f Rl{s)So{s)ds 

Jzkoo 


(5.13) 


ds{Rl{s))ds, 


2 J iboo 

which yields the conclusion. ■ 

Armed with the above lemma, a statement regarding P^ and P^ can now be made. 
Lemma 5.12 When e = 0, 

pi ^ P' pi ^ c' 

^r(l) 2 ^^^ ** 2 ^^' 
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Proof: Since cf)' = 0 when e = 0, a simple observation yields that 

The conclusion now follows from the above lemma. ■ 

Corollary 5.13 For \z\ < M, 

\Pr%-Pr~<p\=0{e), |P+ -P-|=0(C). 


Remark 5.14 By definition, P 14 = and P 34 = - Pfi^ in Lemma \4.3^ 

It is now desirable to compute P^. First, note that 

) = PrUCi,Ct) = 0 , 

and that 


Since 


Prc{Cl,Ct)=RS. 
Prs ~ ~SPrs + + Pj-ui — {eaS — 


(5.14) 


— 2 ea [(6 + aoj) + (di — a)Efi' + (^2 — aa)R^]Pre-, 
upon substitution of the above relations one sees that 

(P±)' = -SPti + 2<1>P± - [eaS - 4.)PP. 

The solution to this equation is given by 

R{z)Pi^g{z) =—ea [ R^{s)S'^{s) ds + [ P(s)<h(s)(P(s)S(s) + 2Pi(s)) ds. (5.15) 

J diOO J drOO 


Using Lemma 5.12 and the fact that P' = RS yields that for bounded z, 

R{z)S{z) + 2P^^{z)=0{e), 

which implies, since <h = 0(e), that the second integral is O(e^). The above argument gives 
the following lemma. 

Lemma 5.15 Let M > 0 be given. Then 

RQ{z)Ppg{z) = —ea f (Po)^(s) ds + O(e^), ze{—oo,M] 

J — 00 

/ OO 

(Po)^(s) ds + O(e^), 2 ;G[—M, 00 ). 

Proof: The conclusion follows immediately from (5.15), taking asymptotic expansions for 
R and S, and using the fact that R' = RS. ■ 

From this lemma one can derive the following three corollaries. 
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Corollary 5.16 Let M > 0 be given. Then for \z\ < M, 

RQ{z)Prs{z) = -A'^ae + O(e^), 

where Prs = Pfg — P^s ■ 


Corollary 5.17 The quantity 9cr“(0) has the asymptotic expansion 

acr“(0) = Nae + 0{e^), 

where N < 0 is given by 


N = 


A'o 


2^o(0) 


Proof: Given the result of Lemma ^.15 , the conclusion follows immediately from the facts 
that 

p-( 0 ) = -S'(0)9,r“(0), 
and that 5'(0) = i?"( 0 )/i?( 0 ). ■ 


Recall the expression given for E'" fd) in Lemma Given the results of Gorollary ^.13 
and Lemma 5.15, a definitive statement can now be made about this quantity. 


Corollary 5.18 Suppose that E"{0) = 0. Then 

L;"'(0) = -{Na + 0{e-^^^))e + O(e^), 

where N > is given by 

N = 3A'25^A2. 

Proof: Substitution of the expression for Prs into the expression for £’^''(0) yields that 

/ CO 

Ro{s)d^Ro{s) ds + 0{e-^^^)e + 0{e^) 

-OO 

= -(3A'25^A20 + 0(e-^4'^))e + 0(e2). 


The fact that the constant N is positive follows immediately from Corollary A.2 


6 Final Arguments 

By Lemma [4.1|, Lemma and Corollary ^.17 it can be seen that 


where 

and 


£"(0) = Cia(6 - C 2 {d 2 - -aa))e‘^ + 0{e^), 

o 


Cl — 2A2A29(jA24 < 0 


( 6 . 1 ) 


c2=-|4^>o, 

dujA24 
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Furthermore, when £^"(0) = 0, by Corollary ^.18 

£'"( 0 ) = -£306 + 0 ( 62 ), 

where £3 > 0 . 


( 6 . 2 ) 


The proof of Theorem is now essentially complete. The result follows immediately 
from the expansions given in equations 0 and (|6.2D , and the fact that £('^^(0) < 0 
(Corollary |2.3| ). The reason the eigenvalues are 0(6) and real follows immediately from the 
fact that £"(0) = 0(6^), while £'"(0) = 0(e). 

The conclusion of Theorem |1.6| follows immediately from the work of Kapitula and 
Maier-Paape |^]. In order to use that work to conclude the existence of multiple pulse 
orbits, all that is necessary is to show that duj(p^(0) 7 ^ 0. This condition is met as a 
consequence of Corollary The fact that the multiple pulse solutions are unstable for 

b < b* follows immediately from the fact that the primary pulse is unstable for b < b* 
(Alexander and Jones Q). The minimal number of unstable eigenvalues also follows from 
that work. 
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A Appendix: Evaluation of Constants 


In the following, 


P = — -auj. 


Proposition A.l The quantities A 2 and A 4 satisfy 

1. A 2 = Ay/Zj-^Tanh~^{\/P) 

VP 

2 . A 4 = - ^A 2 ). 

Proof: Parts 1. and 2. follow immediately from a direct integration, and can be verified 
with the help of Maple V (Release 4). ■ 

Corollary A.2 The functions A 2 and A 4 satisfy 

2 . dujAi = 4:iodojA2. 

Using the fact that 

Tanh“^(x) = - In—iL^ 

^ ^ 2 1-x’ 

Taylor series can be generated for various quantities. 

Corollary A.3 When 0 < /3 < 1 one has the Taylor expansions 


^ fJTl 

1. A 2 = 

2n + 1 


2 . A 4 = ^ 


2n + 3 


3 92 92 . 32 92 .23 

3. A 24 = —{1-—P - 5^/32 - + 

4w^ 3-5'^ 3-52-7^ 3-53-7^ 

. A 8 32 23 ^2 3-1879 ^,^4aa 

4- Ad2 - + 53.72 . 11^ +C’(/^)) 

5. 5^Ad2 = ^(1 + + ^5^77^^ + ^53.72.11^^ + C)(/5"^)) 
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